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ABSTRACT 

We present new VLA C-configuration H I imaging of the Lya-emitting star- 
burst galaxies Tol 1924-416 and IRAS 08339+6517. The effective resolution 
probes neutral gas structures larger than 4.7 kpc in Tol 1924— 416, and larger 
than 8.1 kpc in IRAS 08339+6517. Both systems are revealed to be tidally in- 
teracting: Tol 1924-416 with ESO338-1G04B (6.6' = 72 kpc minimum sepa- 
ration), and IRAS 08339+6517 with 2MASX J08380769+6508579 (2.4' = 56 kpc 
minimum separation). The H I emission is extended in these systems, with tidal 
tails and debris between the target galaxies and their companions. Since Lya 
emission has been detected from both of these primary systems, these obser- 
vations suggest that the geometry of the ISM is one of the factors affecting the 
escape fraction of Lya emission from starburst environments. Furthermore, these 
observations argue for the importance of interactions in triggering massive star 
formation events. 

Subject headings: galaxies: starburst — galaxies: interactions — galaxies: indi- 
vidual (Tol 1924-416) — galaxies: individual (IRAS 08339+6517) 

1. Introduction 

Here, we present new VLA H I imaging of the starburst galaxies Tol 1924—416 [a,6 
(J2000) = 19:27:58.2, -41:34:32; D = 37.5 Mpc (Ostlin, Bergvall, & Roennback 1998)] 
and IRAS 08339+6517 [a,6 (J2000) = 08:38:23.2, 65:07:15; D ~ 80 Mpc (Assuming Hq = 
72 kmsec"^ Mpc^^; Preedman etal. 2001)]. Both galaxies reveal extended tidal structure 
in neutral hydrogen and are clearly interacting with nearby companions. Each appears to 
provide evidence for tidally-induced starburst episodes. 

These two systems were chosen for imaging because they exhibit prominent Lya emis- 
sion (Leitherer etal. 1995; Giavahsco, Koratkar, Sz Calzetti 1996; Kunth etal. 1998; Leitherer 
etal. 2002). Hydrogen Lya is one of the most important diagnostic emission hues in astro- 
physics. It is predicted to be luminous in star-forming galaxies (Chariot & Fall 1993), and 
could potentially be used as an indicator of star formation activity in high-redshift systems. 
However, this application is not straightforward, since the propagation of Lya photons in 
starburst galaxies appears to be a complex process that depends on many factors intrinsic 
to both the starburst itself (e.g., star formation rate, metallicity) and to the surrounding 
galaxy (e.g., geometry, kinematics, and dust content). 
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In a large spectroscopic sample of comparatively high-redshift Lyman break galaxies 
(LBGs; these systems have typical star formation rates > 10 Mq yr~^ and are thus directly 
comparable to more local "starburst" galaxies), Shaplcy ctal. (2003) find various correla- 
tions between the emergent Lya profile (i.e., strengths of absorption or emission) and char- 
acteristics of the ISM. In particular, dustier galaxies have less prominent Lya emission, or 
equivalently, redder UV continuum slopes. Lya emission strength is shown to anti-correlate 
with the kinematic offset implied by the redshift difference between Lya emission and low- 
ionization interstellar absorption lines; that is, systems with larger outflow speeds usually 
show less prominent Lya emission. Finally, systems with stronger Lya emission show lower 
star formation rates, on average, than systems with weaker Lya emission. These correla- 
tions suggest that many factors can contribute to the escape fraction of Lya from starburst 
regions. 

One of the most important factors that influences Lya propagation is resonant scatter- 
ing. It is well known that due to high resonant scattering by neutral hydrogen in the ISM, 
Lya photons can be attenuated by even small amounts of dust. It is then expected that 
only young, relatively dust-free galaxies should be prodigious sites of Lya production. Low- 
metallicity starburst galaxies in the local universe may be considered nearby analogs to such 
objects which are expected in greater numbers at higher redshifts. Thus it was surprising 
that HST observations of the most metal-poor galaxy known, I Zw 18, showed only damped 
Lya absorption and no emission (Kunth etal. 1994). In stark contrast, the more metal- rich 
starburst galaxy Haro2 showed prominent Lya emission (Lequeux etal. 1995). 

The effects of resonant scattering can be lessened if there is a velocity difference between 
the site of Lya emission and the surrounding material from which resonant scattering can 
occur. This will Doppler shift the Lya photon out of resonance and increase its chance 
of escaping the surrounding ISM. While this process seems to be complicated (note, for 
example, that higher outflow speeds suggest less prominent Lya emission profiles in the 
sample of Shaplcy ct al. 2003), it appears that an important factor is the presence or absence 
of an outflow (see below). Indeed, Shapley etal. (2003) measure an average Lya offset > 
300kms~^ for the sample of Lya-emitting LBGs. These observations and the studies of 
Giavalisco etal. (1996) and Kunth etal. (1998) argue for the importance of ISM geometry 
and kinematics in affecting the Lya escape fraction. 

The studies of Kunth etal. (1998), Tenorio-Tagle etal. (1999), and Mas-Hesse etal. 
(2003) have examined the characteristics of the ISM kinematics which appear to influence 
the propagation of Lya photons in star-forming galaxies. If static, homogeneous neutral gas 
with column densities > 10^^ cm~^ shields the ionized gas, no emission will be detected. The 
resonant scattering of the Lya photons will lead to increased probability of destruction by 
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any dust which is present. On the other hand, there may be diffuse Lya emission which 
is detectable on sightlines not coincident with the sources of UV photons. Similarly, if the 
areal coverage of the neutral gas is not uniform but clumpy, some Lya emission may be 
detectable on favorable sightlines. Finally, if the velocity structure of the neutral gas is not 
static but rather outflowing from the ionizing regions (outflow velocities > 200kms~^), Lya 
photons to the red of 1216 A can escape and Lya emission may be significant. This explains 
the strong Lya emission detected in some starburst galaxies with complete spatial coverage 
by neutral gas which is also comparatively rich in both metals and dust. 

These observations were undertaken with these arguments in mind. Our group has 
obtained new HST/ACS imaging of the Lya line in a small sample of galaxies (sec Kunth 
et al. 2003), two of which are Tol 1924-416 and IRAS 08339+6517. While these data are only 
a small segment of the sample that was imaged in Lya, they will aide in the interpretation 
of the Lya data. We note that our H I imaging sample does not include systems showing 
only Lya absorption, so we must await a larger sample to draw definitive conclusions on 
the correlation between H I kinematics and the appearance of Lya emission. However, the 
results presented here support the scenario that the geometry and kinematics of the ISM are 
important (although not the only) factors that affect the escape probability of Lya photons 
from massive star formation regions. 



2. Observations and Data Reduction 

The NRAO Very Large Array^ (VLA) was used to obtain H I spectral line data for 
these two systems, as part of observing program AC 654. Tol 1924— 416 was observed in 
the BnC configuration for 180 minutes on 2002 September 30, while IRAS 08339+6517 was 
observed in the C configuration for 370 minutes on 2002 December 8 & 9. The correlator 
configuration yields a bandwidth of 3.125 MHz, with 64 channels separated by 48.8 kHz 
(10.5 kms~^). The data were reduced and calibrated using standard methods in the AIPS 
environment. The data were imaged using uniform weighting; the single-plane rms noise is 
1.3 mJyBeam-^ for Tol 1924-416 and 0.40 mJyBeam-^ for IRAS 08339+6517. With the 
final beam sizes (26"xl6" for Tol 1924-416, and 21"xl6" for IRAS 08339+6517), this noise 
level implies single-plane 3a detection thresholds of 2.3x10^° cm~^ and 6.9x10^^ cm~^ for 
Tol 1924-416 and IRAS 08339+6517, respectively. 



^The National Radio Astronomy Observatory is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, Inc. 
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3. Hi Data: Extended Tidal Structure 

With these data we detect extended tidal structure in each of these systems. In both 
cases H I gas is detected in the target system, in a neighboring galaxy, and in tidal material. 

3.1. Tol 1924-416 

In Figure 1 we present contours of the zeroth-moment (total column density) image 
overlaid on a DSS optical image of Tol 1924— 416, as well as an intensity-weighted velocity 
field. We detect a total integrated fiux of 12.6±1.58 Jykms"^ from the Tol 1924— 416 sys- 
tem. This compares well with the single-dish value of 11.4 Jykms"^ derived from Parkes 
observations in the HIPASS Survey^ (Barnes etal. 2001), suggesting that we have missed 
little or no flux due to the lack of zero-spacing interferometric elements. At the adopted 
distance of 37.5 Mpc, this corresponds to a total H I mass of (4.2±0.5) xlO^ M©. 

Inspecting the H I data cube allows us to measure the masses of neutral gas clearly 
associated with both Tol 1924—416 and with the companion ESO 338— IG04B. In the former 
we detect a total mass of (1.4±0.2)xl09 Mq, and in the latter we find (9.3±1.2) xlO^ Mq. 
The remainder of the system mass, (1.9±0.2) xlO^ Mq or ~ 40%, is tidal material (see 
Figure 1). 

The velocity structures of Tol 1924-416 and ESO 338-IG04B have been studied in the 
Ho; fine by Ostlin etal. (1999) and Ostfin etal. (2001). The agreement between the optical 
and radio velocity fields is very good. The radio data confirms the weak large-scale ordered 
rotation in Tol 1924—416, and the very ordered rotation in the companion system. Note that 
the H I associated with the Lya-emitting system (Tol 1924—416) is spread out over ~ 80 
km s~^ in radial velocity, suggesting that there is not a static screen of neutral gas along the 
line of sight to the starburst region. 

Tol 1924— 416 shows a variety of interesting characteristics in optical imaging. Ostlin 
et al. (2003) used multicolor HST images to derive the cluster formation history. The elevated 
star formation rate (~ few MQvr^-*^) has been present for 40 Myr, and shows evidence for 
propagation throughout the disk. There appears to be little to no reddening in this system, 
in agreement with expectations based on the sample of LBGs in Shapley et al. (2003). While 
each of these parameters likely influence the propagation of Lya photons through the ISM, 



^The Parkes telescope is part of the Australia Telescope which is funded by the Com- 
monwealth of Australia for operation as a National Facility managed by CSIRO; see 
http: / / www.atnf.csiro.au/research/multibeam/ 
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our observations suggest that the ISM kinematics are also an important factor in this process. 

The angular distance between the two galaxies is ~ 6.6', corresponding to a minimum 
deprojected distance of ~ 72 kpc. The total radial velocity extent is ~ 150 kms~^ {^sys — 
2830 kms~^), but components exist in both systems at the same velocities. Hence it seems 
hkely that these two systems may now be gravitationally bound. Regardless, the extended 
nature of the H I suggests that the present starburst in Tol 1924—416 likely was triggered by a 
strong gravitational interaction with ESO 338— IG04B, and that the different star formation 
rates in these two otherwise similar systems [comparable H I masses (see above) and stellar 
masses (see Ostlin et al. 2001)] are attributable to the recent tidal interaction between them. 



3.2. IRAS 08339+6517 

In Figure 2 we present contours of the zeroth-moment (total column density) image over- 
laid on a DSS optical image of IRAS 08339+6517, as well as an intensity-weighted velocity 
field. We detect a total integrated fiux of 3.68±0.46 Jykms"^ from the IRAS 08339+6517 
system; this compares well with the single-dish value derived from Nangay observations, 
3.81 Jykms^^ (Martin etal. 1991), again suggesting negligible flux loss due to lack of zero- 
spacing basehnes. At the adopted distance of 80 Mpc, this corresponds to a total H I mass 
of (5.6+0.7) xlO'^ Mq. 

We flnd (1.1+0.2) xlO^ M© of H I to be associated with IRAS 08339+6517. Its compan- 
ion system, 2MASX J08380769+6508579, is found to have a total H l content of (7.0+0.9) x 10^ 
Mq, leaving the remaining (3.8+0.5) x 10^ Mq to reside in tidal material between the two 
systems. This implies that the bulk of the neutral gas (~ 70%) in this interacting system 
has been removed from one or both galaxies. Interestingly, the peak single-plane column 
density is found in tidal material and is not associated with either galaxy. 

Mas-Hesse etal. (2003) have shown from the analysis of HST/STIS high-resolution UV 
spectroscopy that in front of the central ~ 5 kpc (13") of this galaxy a column density ~ 
-|^Qi9.9 neutral gas is outflowing with a velocity of around —300 kms"-*^. This flow of 

material leaving the central area of IRAS 08339+6517 could be feeding the region between 
both galaxies. It is interesting to note that the neutral gas in this intermediate region is 
approaching us with an average velocity of around —100 kms~^. However, no neutral gas 
is detected in the STIS observations at the systemic velocity of the galaxy {Vsys = 5750 
kms~^) derived from these radio observations. This implies the existence of flner details in 
the velocity structure than can be inferred from these H I observations (e.g., an expanding 
shell powered by the starburst; see Mas- Hesse etal. 2003). 
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The distance of IRAS 08339+6517 (80 Mpc) precludes detailed optical studies of the 
stellar or cluster populations. Inferences about other intrinsic properties of the galaxy must 
therefore rely on spectral data. Gonzalez Delgado etal. (1998) measure an Ha flux corre- 
sponding to a star formation rate of ~ 8 Mq yr~^ (applying the prescription of Kennicutt, 
Tamblyn, & Congdon 1994). The measured reddening reaches values as large as E(B— V) = 
0.5 mag., suggesting substantial amounts of dust in this system. Again we note that these 
factors likely contribute to the propagation of Lya photons. However, since this starburst 
region appears to suffer from substantial reddening, the importance of ISM kinematics in 
Doppler shifting the Lya photons out of resonance appears to be pronounced. 

For a measured angular separation of 2.4', the implied minimum separation is 56 kpc. 
The H I radial velocity extent is ~ 300 kms~^; this suggests that the neutral gas surrounding 
IRAS 08339+6517 is highly turbulent, having been powered by the release of mechanical 
energy in the central starburst. At a smaller spatial scale the turbulence is still higher, as 
derived from UV spectroscopy (see Mas- Hesse etal. 2003). 

While we can not definitively state that these galaxies now form a bound system, the 
effects of the gravitational interaction are pronounced in IRAS 08339+6517. Interestingly, 
this starburst system shows two peaks of Lya emission in our ACS imaging, and may be the 
remnant of an earlier merger in its own right. On the other hand, the comparatively under- 
luminous nature of 2MASX J08380769+6508579 (i.e., no obvious massive starburst episode 
in the current epoch) suggests that the effects of the interaction have not been as severe for 
the secondary system. Again, the less-massive system appears to retain a strong component 
of solid-body rotation. IRAS 08339+6517 also appears to be undergoing organized rotation, 
but the neutral gas is spread out along the line of sight over ~ 80 kms~^. 



4. Conclusions 

VLA H I imaging of the starburst galaxies Tol 1924— 416 and IRAS 08339+6517 has 
been presented. These two systems are remarkably similar in H I content, mass, and current 
evolutionary state. In each, we find extended neutral gas between the target and nearby 
neighbors, suggesting that interactions have played an important role in triggering the mas- 
sive starbursts in the primary galaxies. The close proximity of the companions suggests 
that the interactions were recent, and the similar velocities of both primary and secondary 
galaxies argues that these systems may end up gravitationally bound. 

Since both primary systems are intense Lya emitters, these data support the inter- 
pretation that the ISM kinematics are an important mechanism that affects the escape of 
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Lya photons from starburst regions. Combined with the results on LBGs (Shapley et al. 
2003), these data have immediate imphcations for the use of the strength of Lya emission 
in determining star formation rates, since the results will be dependent on the geometry of 
the ISM and other factors (e.g., dust, star formation rate) and not on properties inherent 
to the starburst being considered. Further H I observations of Lya-emitting galaxies (and, 
conversely, of starburst systems with no apparent Lya emission) are certainly warranted to 
further explore the role of the ISM in regulating the escape of Lya photons from starburst 
environments. 
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Fig. 1. — (a) DSS image of Toll924— 416, overlaid witli contours of tlie H I zerotli-moment 
image. Contours correspond to column densities of (7.3, 29, 51, 73)xl0^° cm~^. Each 
galaxy is labeled; beam size is shown at bottom left, (b) Intensity-weighted velocity field of 
Tol 1924—416. From this figure it is apparent that H I is being removed from one or both 
systems. Also, there remains a component of solid-body rotation within the optical extent 
of both galaxies. Beam size is labeled at lower left. 
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Right Ascension (J2000) 

Fig. 2. — (a) DSS image of IRAS 08339+6517, overlaid with contours of the H I zeroth- 
moment image. Contours correspond to column densities of (5.5, 15, 24, 33, 42, 51)xl0^° 
cm~^. Each galaxy is labeled; beam size is shown at bottom right, (b) Intensity- weighted 
velocity field of IRAS 08339+6517. From this figure it is apparent that H I is being removed 
from one or both systems. The companion galaxy appears to retain a component of solid- 
body rotation in neutral gas; clear signs of rotation are less prominent in IRAS 08339+6517, 
however, suggesting that this interaction has completely disrupted the neutral gas in this 
system. Beam size is labeled at lower right. 



